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Abstract 

The Wilson-Dirac matrix for 577(2) with I = \ fermions is written in an explicitly 
real form. The basis change relating it to the conventional form in which the matrix has 
complex entries is also given. Some applications are presented. 
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The Wilson-Dirac operator for SU(2) with fermions in the fundamental representation 
on a lattice of size L 4 is usually represented by an 8L 4 x 8L 4 complex matrix. Below, I 
first present the Wilson-Dirac operator as a real matrix of the same size. Next, I write 
down the basis change relating it to the complex form. I end with some comments. 

To construct the Wilson-Dirac operator we need, at each lattice site, real linear com- 
binations of elements of an algebra generated by 1, {7^, /j, = 1, 2, 3, 4}, {icrj, j = 1, 2, 3}, 
acting on an 8 x 8 dimensional space. 1 is the 8x8 unit matrix. The 7^ and aj (no 
i-factor) are hermitian. The following relations define the needed seven matrices up to a 
choice of basis: 



hvivj] = , iajicTk = -e jk iiai , (iaj) 2 = -1 , {7^7^} = 25^. 



(1) 



The real form of the Wilson-Dirac operator is now defined by exhibiting seven real 
8x8 matrices obeying (1). To distinguish these matrices from the more standard, complex, 
versions I denote them by the corresponding capital letters. For practical purposes I looked 
for a form in which the T M matrices are sparse, have simple numerical entries, and have 
an off-diagonal structure characteristic of a chiral basis. 

The T M have the following structure in terms of 4 x 4 blocks: 
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where denotes the transpose of p M . The iT,j have the following structure in 4 x 4 blocks: 

s] = -1. (3) 

The explicit 4x4 matrices we need are listed below: 
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The above representation is related by a basis change to the usual, complex, version. 
More precisely, the 7^ matrices are taken in a chiral basis as complex 4x4 matrices: 

Ik = 01 ® a k , fc = l,2,3, 74 = 02 <S> 1, (5) 

where the 2x2 matrices at are the standard Pauli matrices acting on two dimensional 
complex space. The Pauli matrices are also chosen for the action on the group index. I 
replace the single four valued spinorial index in the usual representation by two two- valued 
indices matching the direct products of equation (5). Thus, a fermion field in the usual 
representation carries four indices: The first three indices take two values and label the 
eight dimensional space per site. The last index takes L 4 values and labels the site. The 
first pair of indices are Dirac and the third index is SU (2)-color. 
Acting on these fermion fields we define the unitary matrix $: 

§ = ]-u®[a® (3 + [3 ® a + [3® (3 - a® u]®l, (6) 
where a and (3 are given by 

-c !)■ "-a d- ^ 

The four 7^ matrices turn into the T M matrices under conjugation. The latter are writ- 
ten below using three two valued indices rather than the one eight valued index employed 
equations (3,4): 

$(7l ®1® l)$ f = Ti <g) 1 = -6(g) £ <g) 1 <g) 1 
$(72 <8> 1 <S> l)$ f = T 2 ® 1 = -£ <S> (7i ® £ ® 1 

(8) 

$(73 ® 1 ® 1)$ T = T 3 <g> 1 = £ <8) 03 ® £ <S> 1 

$(74 ® 1 <s> i)$ f = r 4 <8) 1 = -ai <g> 1 <g> 1 <g) 1, 

where 

e = -^0 2 = (J "q 1 ), e 2 = -l. (9) 

Similarly, the iaj matrices become under conjugation the iHj matrices of equations 
(3,4): 

$(1 <g> 1 <g> (ioi) <g> l)$ f = iEi <g> 1 = -1 <8) 1 <g) £ <8) 1 

$(1 <g> 1 <g> (icr 2 ) <S> l)$ f = iE 2 <g> 1 = -1 <g> e <g> <ri ® 1 ( 10 ) 

$(1 ® 1 <g> (zcr 3 ) ® l)$ t = iE 3 (8) 1 = 1 <g) £ ® a 3 ® 1. 

The reality of the Wilson-Dirac operator should be useful in simulations of SU(2) 
gauge theories. From eq. (4) we see that the action of the matrices amounts to permu- 
tations plus sign switches. The projector combinations one encounters in the Wilson-Dirac 
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operator, 1±T M , are easily efficiently implemented because the matrices pi.,2,3 are antisym- 
metric and P4 is symmetric. When compared to QCD one would expect the computational 
gain in fermion manipulations (for example in the inversion of the Wilson-Dirac matrix) 
when restricting the gauge group to 577(2), but staying with a complex representation, to 
be by a factor of about |. I have implemented the above real form in a parallelized high 
level Fortran 90 code and a comparison has given a gain of about 6. A beta version of the 
package I used is available upon request. A real form of the Wilson-Dirac operator can be 
obtained for any non-integral isospin, but not for integral isospin. 

The above basis is particularly useful for the overlap [1]. The matrix 75 is diagonal 
before and after the $ transformation. The Hamiltonian matrix used in the overlap, H, 
is given by T 5 Dw, where Dw is the Wilson-Dirac matrix. H is both real and hermitian, 
hence symmetric. All the eigenvectors of H can be chosen real. Thus, the overlap itself is 
real. This means that a theory representing a single Weyl doublet in interaction with the 
577(2) gauge field, has a real fermion determinant on the lattice just like in the continuum. 
Under gauge transformations sign switches can occur. This is how Witten's anomaly [2] 
will be reflected on the lattice. On the basis of [3] these sign changes can be associated 
with degeneracies of H of codimension 2 in this case (conical degeneracies [4]). Such 
degeneracies are known to occur in the continuum by Witten's argument [2]. 

What we need here is made explicit in what follows: Consider first the massless Dirac 
operator, in Euclidean continuous space, made hermitian and real. We take space as a 
torus of size / in each direction. Let us label the space-time points by four vectors £. Fixing 
£4, we have a three torus which can be mapped into 577(2) with unit winding [5]: 

ina>42,43j= -, l-L-U 

+ EUli " ^(x^)]} 2 + E' =1 [sm(^^)] 2 

This choice preserves the symmetry under discrete rotations compatible with the three 
torus. Now, following a formula due to Goldstone, presented in [6], we introduce a gauge 
transformation depending on all four coordinates: 

g(0 = e < ^Gt(e 1 ,6,6)e- < ^G(ei,6,e 3 ). (12) 

This map cannot be deformed to unity. [2] shows that as the vector potential associated 
with Q is deformed to zero a conical singularity will occur in the massless Dirac operator 
somewhere along the way. 

It is trivial to put Q on the lattice, and to again define a link by link interpolation 
to unit link variables everywhere. Now, in addition to the deformation parameter, we 
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also have the mass parameter m in D\y (alternatively known as the hopping parameter 
k) at our disposal. To be definite, we choose m = — 1 (corresponding to k = |) and start 
increasing it towards zero. For positive m (corresponding to k < |), H has a gap at zero 
[1] and nothing can happen. However, the smoothness of the gauge configuration in (13) 
implies that in the limit of infinite L the continuum conical singularity must appear at 
zero m. It cannot come from nowhere, so it must be somewhere in the plane spanned by 
the real deformation parameter and negative m. 

For SU(2), the reality of the Wilson-Dirac operator also extends to truncations of the 
overlap, which are almost massless lattice versions of continuum vector-like gauge theories 
[7]. In the context of chiral gauge theories we recall that while SU(2) gauge theory with a 
single Weyl fermion of isospin ~ is inconsistent because of the global anomaly, if the isospin 
is |, the problem disappears [2]. Taking one Weyl fermion of isospin | is probably the 
simplest consistent four dimensional chiral model with a real fermion determinant. The 
reality is similar to the 11112 model in two dimensions [8] which yielded a successful test 
of the overlap. On the other hand, the four dimensional model has not been exactly solved 
in the continuum and its properties are unknown. 
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